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Abstract. The use of Reduced Proper Motion in identifying isolated white dwarfs 
has long been used as a proxy for the absolute magnitude in a population with known 
kinematics. This, however, introduces a proper motion detection limit on top of the 
existing photometric limit. How the survey volume is hampered by this extra parameter 
is discussed in Hambly et al. (2012). In this work, we discuss some robust outlier 
rejection methods in order to minimise the proper motion limit and hence maximise the 
survey volume. The generalised volume, corrected for the distance of the Sun from the 
Galactic Plane, is integrated explicitly. 


1. Introduction 

Main sequence (MS) stars with masses of less than 8 M© will end up as white dwarfs 
(WDs). This mass range encompasses the vast majority of stars in the Galaxy. Thus 
a WD, which is a degenerate core left behind from its progenitor, is the most com¬ 
mon end-point of stellar evolution. Nuclear burning is negligible at this stage, so WDs 
cannot replenish the energy they radiate away. Hence, the luminosity and temperature 
decrease monotonically with time. The WD luminosity function (LF) was first used 
as a cosmochronometer half a century ago. Given a finite age of the Galaxy, there is 
a temperature beyond which the oldest WDs have not had time to reach which would 
translate to a sudden downturn in the WDLF. Further to cosmochronometry, a WDLF 
can be inverted to provide the star formation history (Rowell 2013). 

2. Pan-STARRS 

Pan-STARRS-1 (PSl) is a 1.8 m optical wide-field imager developed by the University 
of Hawaii (See Metcalfe et al. (2013) and references therein for details). The PSl 37r 
Steradian Survey covers the sky north of declination -30° with five broadband filfers 
designafed as gpi, rpi, ipi, Zpi and ypi spanning 400-1000 nm. The sky has been imaged 
60 limes on average in Ihe 4-year survey, allowing PS 1 fo delermine accurate proper 
molions wilh small epoch differences. The lesl dala covers Ihe sky belween declination 
0° and 7.5°. The selection criteria are: 

1. A minimum of 1.5 years epoch difference 

2. Proper molion wilh lOcr confidence 

3. Deleclion in all five filters 
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4. Good morphology flags 

2.1. Proper Motions 

In the PSl reduction pipeline, the proper motions and the parallaxes are calculated 
simultaneously for all objects having sufficient coverage in parallax factor. However, 
most objects would not have detectable parallaxes such that some objects can have 
either spurious proper motions or unrepresentative uncertainties and values for the 
proper motions. Furthermore, in a proper motion limited sample, the sampling volume 
scales between yp- and yp depending on the population (Hambly et al. 2012). Therefore, 
a smaller lower proper motion limit would greatly increase the survey volume. We 
have investigated some robust algorithms by rejecting outliers to improve the quality of 
proper motion, hence to reduce the proper motion limit to be adopted: 

1. Fitting by minimising absolute deviation 

With this method, data points are unit weighted. It is found to be ineffective at 
handling clustered data. 

2. Jackknife method 

This would be the ideal option, however, there is not enough computing power to 
go through this process (jackknifing 10 points out of the 60 epoch measurements 
for 10^ objects would require 10^^ calculations). 

3. Iterative outlier rejection 

This method identifies data points lying outside 3cr from the best fit solution as 
outliers. A new weighted-least square solution will be found and this process 
continues until no more data points are rejected. This method is much faster than 
jackknife, but it is more sensitive to data further away from the centroid of the 
solution. 

4. Improved iterative outlier rejection 

The improved version has the tolerance level based on the propagation of errors 
which accepts larger deviations when the data points are further away from the 
centroid. 

The best fit line can be described by the function 


f = X+flt 


( 1 ) 


and the uncertainty is given by the standard propagation of errors using the covariance 


matrix 




2 , 2 ^" 


( 2 ) 
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3. Methods 

3.1. Model Atmosphere 

The synthetic colour of both DA, DB and mixed hydrogen-helium model atmospheres 
in the Pan-STARRS colours are provided by Dr. Bergeron (Holberg & Bergeron 2006; 
Kowalski & Saumon 2006; Tremblay, Bergeron & Gianninas 2011; Bergeron et al. 
2011) which were based on the most recent calibrations. At this early stage of analysis, 
it is assumed that all WDs have pure hydrogen atmosphere and have surface gravity 
log(g) = 8.0 in order to fit the effective temperatures and the distances simultaneously. 

3.2. Reduced Proper Motion 

Since WDs have much smaller radii than MS stars at the same temperature, the WD 
cooling sequence is a few magnitudes fainter than the main sequence. An efficient way 
to select WD candidates is to use reduced proper motion (H) (Kilic et al. 2006; Harris 
et al. 2006). By using proper motions as proxy-parallaxes, one can obtain the reduced 
proper motion, which is analogous to the absolute magnitude: 


Hrpi ^5 + 5 logQi) + rpi 


(3) 


Using proper motion and the photometric parallax, the tangential velocities, Vtan, can be 
deduced. This is an important quantity in deriving the distance limit due to the proper 
motion limit. By holding Vtan constant, one can determine how far the object can be 
placed before its proper motion would drop below the lower proper motion limit. 

3.3. Generalised Schmidt Estimator 

In the TT^— method, the contribution of each object to the LF is weighted by the inverse 

» max 

of the maximum volume in which the object could be observed by the survey. However, 
this technique assumes objects are uniformly distributed. In reality, stars in the solar 
neighbourhood are concentrated in the plane of the disk so the effects of space-density 
gradient have to be corrected. This led to the development of the generalised volume 
Vgen (Stobie, Ishida & Peacock 1989) which is calculated by integrating the appropriate 


stellar density profile — = exp(-^) = exp(-^^^^^) along the line of sight. 



(4) 


where ;\^(vtan) is the discovery fraction of the sample from the lower tangential velocity 
limit, D is the size of the solid angle of the survey, dmin and d^ax are the distances 
limits set by the bright and faint detection limits as well as the high and low proper 
motion limits, H is the scale height of the Galactic disk profile. By further taking into 
account of the distance of the Sun from the Galactic Plane, the density profile becomes 
^ = exp and the integral can be solved analytically: 



V H 

dmax sm(b)+Zo 

H 


dmax sin(b)+Ze 


(2^3 + 2dmax^^ + d^ax^) > if [dmax sin(b) -I- Zq] > 0 
) (2^3 - 2dmax^^ + d^x^ > Otherwise, 


'max' 


max 


where ^ = (H/ sin (b)). 


(5) 
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4. Future Work 

The final data release of PSl is scheduled to be in mid 2015. The data quality is ex¬ 
pected to increase with the improvement in the reduction pipeline and the increase in the 
maximum epoch difference. The total number of WD candidates with Vtan > 40 km s“' 
is expected to be about 40000 (Hambly et al. 2012). With sufficient objects, it is possible 
to untangle the thin disk, thick disk and stellar halo which did not anive at statistically 
confident results in previous work (Rowell & Hambly 2011). Furthermore, the signif¬ 
icant increase in the survey volume allows the inversion of the luminosity function to 
recover the star formation history of the Galaxy with greater precision. 
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